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We study the zero-temperature phase diagram and fractional excitations when a thin film of 3D 
topological insulator has two competing masses: T- symmetric exciton condensation and T- breaking 
Zeeman effect. Two topologically distinct phases are identified: in one, the quasiparticles can be 
viewed as in a quantum spin Hall phase, and in the other a quantum anomalous Hall phase. The 
vortices of the exciton order parameter can carry fractional charge and statistics of electrons in both 
phases. When the system undergoes the quantum phase transition between these two phases, the 
charges, statistics and the number of fermionic zero mode of the excitonic vortices are also changed. 
We derive the effective field theory for vortices and external gauge field and present an explicit 
wave function for the fermionic zero mode localized at the excitonic vortices with or without orbital 
magnetic field. The quantum phase transition can be measured by optical Faraday or Kerr effect 
experiments, and in closing we discuss the conditions required to create the excitonic condensate. 

PACS numbers: 71.35.-y, 71.10.Pm, 73.20.-r 



I. INTRODUCTION 

The surface states of three-dimensional topological in- 
sulators (TI) are two-dimensional metals with an odd 
number of Dirac points enclosed by the Fermi surfacd^^'^. 
They avoid the so-called fermionic doubling problem and 
are known in several case^SHU iq support topological cor- 
related states when different types of interactions are 
added. The example motivating this paper is the topo- 
logical exciton condensate formed by the Coulomb inter- 
action between top and bottom surfaces of a thin film 
of topological insnlatoi'^. The exciton condensate on the 
TI surfaces carries electronic charge i^e at its vortices, 
i.e., half of the charge of the underlying electrons. This 
thin TI exciton condensate should be compared to that 
in a graphene bilayer, which has similar relativistic elec- 
tronic structure to a TI surface. The fractional prop- 
erties are washed out by the spin/pseudospin degenera- 
cies in graphene. For example, the Kckulc vortex on 
graphene can support fractional quantum number^^ltiSl^ 
but the fractional excitation should be multiplied by 4 
(or 2, depending on how the degeneracies are lifted) and 
results in integral quantum numbers as for electrons. (See 
also the charge fractionalization on the kagome lattice^^ 
at filling 1/3) 

In the thin film, we have at least two Dirac fermions 
overall, one from the top surface and one from the bot- 
tom, and hence the system looks "ordinary" in the sense 
that we have an even number of fermionic degrees of free- 
dom. However, it still has only half as many degrees of 
freedom as graphene, and it is known from previous work 
that the Topological Exciton Condensate(TEC) is dis- 
tinct from an ordinary BEC of excitons. It is fractional 
in the sense that it supports charge zt^e at its vortices^. 
This fractional quantum number of TEC results from the 
topological structure of the TI. Now the TEC is T- sym- 
metric but generates a mass at the Dirac points. One way 



to motivate the present work is to compare this effect to 
the easiest way to open up a gap at the Dirac points of 
TI surface, namely, to break T-symmetry on the surface. 
There are a few natural ways to break T- symmetry and 
open up gaps at Dirac points: coating magnetic mate- 
rials, or imposing magnetic fields etc. Thus if we intro- 
duce T- symmetry breaking perturbation and condense 
T- symmetric excitons on TI surfaces at the same time, 
then we expect that there should be an interesting in- 
terplay between TEC mass and T- symmetry breaking 
masses. We find a quantum phase transition (QPT) be- 
tween two phases that can be labeled by the behavior of 
their quasiparticles: Quantum Anomalous Hall (QAH) 
and Quantum Spin Hall (QSH) phases. 

Here, we will first consider a thin TI film coated by fer- 
romagnets (FM) on both sides (Figjl]). Then the system 
is gapped due to the T-breaking perturbation of the orig- 
inally gapless surface. The electrons on the TI surface in- 
teract with the FM magnetization via Zeeman coupling 
Jhc ■ M. The ^-component magnetization has a spe- 
cial role: it opens up a gap at the Dirac point with mass 
Jh\Mz \ sgn{Mz). This coupling is different from the case 
of graphene, where the pseudospins do not interact with 
magnetization and the Zeeman effect does not open up a 
gap. Other components of M act as external gauge fields 
coupled to the Dirac fermions. When the system is fully 
gapped, we can safely integrate out fermions and obtain 
effective field theories for the gauge fields. We will de- 
velop the field theories for the fluctuating gauges and vor- 
tices in the two distinct phases. Even though the physics 
is most transparent when only the Zeeman coupling is 
considered, the same physics applies even when we con- 
sider an external z-directional magnetic field. The key 
difference in this case is that we have Landau levels (LL) 
instead of a Dirac spectrum. Despite this, we still have 
fractional vortices and a quantum phase transition when 
the Zeeman energy is comparable to TEC mass. 
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FIG. 1: A. Illustration of a thin topological insulator film 
sandwiched between ferromagnetic layers. To see the quan- 
tum phase transition between quantum spin Hall (QSH) and 
quantum anomalous Hall (QAH) phases, we need to align the 
magnetism along the same direction on both layers. B. Alter- 
nate experimental scheme. We sandwich a thin ferromagnetic 
layer between two bulk topological insulators. If we put an 
insulator with small dielectric constant e instead of ferromag- 
nets, we expect to observe an increased critical temperature 
for the superfluid - insulator transition into the excitonic con- 
densate. 



ferromagnet. We have two Dirac cones indexed by a = 
1, 2 (layer index) and short-ranged Coulomb interactions 
U,V. Here we specialize to the case ^ My <^ M^- 
We include external EM probe fields Ai which couple to 
the Dirac fermions and ignore the tunneling between lay- 
ers, as this decays exponentially in thickness whereas the 
Coulomb energy decays only algebraically. For simplicity 
we set the chemical potential /i = (however, the topo- 
logical properties studied here should be valid whenever 
the chemical potential is in the gap) and start from the 
weak Coulomb interaction where no phase transition is 
expected: 
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Here a are the Pauli matrices for (physical) spin. Ignor- 
ing U = (where I is the system-dependent length of 

short-ranged Coulomb interaction) and V = ^ (where d 
is the width of the thin film), we can easily identify the 
eigenstates and eigenvalues of H. This is reasonable as 
U and V will only renormalize the Dirac velocity when 
they are small. From here onwards, we rescale uj? = 1 as 
this should not change topological properties. With this 
in mind, we rewrite the Hamiltonian in a more familiar 
format H = F^HF where F is the four-spinor for the 
system: 



The following section studies the electronic band struc- 
ture and shows that there is a QPT when the Zeeman 
mass is comparable to the exciton mass. At the QPT, 
the gap vanishes and the topological invariant (Chern 
number) of the bands is changed, which is why we re- 
fer to the transition as between QSH and QAH phases. 
In Section III we obtain an effective field theory that 
captures the topological properties. From this we can 
derive information about the change across the transi- 
tion in the charge, statistics and the number of bound 
zero modes in a vortex. We next proceed to the uniform 
magnetic field case and find generally similar physics as 
in the Zeeman case. We also present the solution of the 
fermion zero mode at the vortices. In QSH, there is only 
one fermionic zero mode which disappears at the quan- 
tum phase transition to QAH. However, it's difficult to 
obtain the exact solution for the fermionic zero mode in 
the uniform orbital magnetic field but we show that we 
can manage to obtain analytic solutions for some limiting 
cases. In closing we discuss experimental detection of the 
phase transition and conditions required to observe the 
exciton condensate. 
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There are two important points to note about this Hamil- 
tonian H. First, the Hamiltonian has the spectrum with 
the gap 2\JhMz — m\ at the Dirac point. 
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The spectrum is doubly degenerate due to the layer in- 
dex. Moreover, from the matrix form([2| of H, we identify 
both up and down layers as in a QAH phase. To see this, 
let us take only the upper half of the block Hamiltonian 
H„ such that 
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Now, we compute the winding number for H„ = p • cr + 
TOtT^. We identify the vector V„ = {px,Py,ni) from H„, 
and compute the topologically invariant winding number 
(the first Chern Number) for V by evaluating 



II. MICROSCOPIC MODEL AND BAND 
STRUCTURE 

We start from a microscopic model of a thin TI with 
a T-breaking Zeeman mass induced by proximity to a 
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where V = V/|V|. The winding numbers for H„ and 
are sgn{m)^ which reflects the origin of the anomalous 
quantum Hall effect for the Dirac fermions. Taken as a 



whole, the system is in a QAH regime where the Hall 
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conductance is quantized as ^. 

So far, we have neglected the Coulomb interactions as 
they only renormalize parameters when they are small. 
However, when the Coulomb interaction increases, there 
can be phase transitions to a stripe phas^^^l for strong 
U , or exciton condensate^ for strong V . We are inter- 
ested here in the TEC which is T- symmetric but still 
opens up a gap. To obtain TEC, we treat V in terms 
of Mean Field Theory(MFT) in the exciton order pa- 
rameter and ignore U . We do the MET for Vnin2 = 
f^j^apj2 ^ ^^^Yi M^P =< /i/2t > (where and 

/^p are two-component Dirac fermions for layer 1 of spin 
^-component a and 2 of spin z-component p). Here, dif- 
ferent M describes different exciton condensates but we 
will choose an M which opens up a gap at the Dirac 
points and is T- symmetric^. The order is described by 
M = Ar^cr" where r are Pauli matrices acting on layer 
index, and a acts on spin as before. With the exciton 
order parameter, we have the following second-quantized 
mean-field Hamiltonian. 
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In this form, the QPT between QAH and QSH phases is 
not obvious. To clarify the physics, we use the matrix 



forms of this Hamiltonian H 
QAH case. 
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The first observation is that the spectrum of the Hamilto- 
nian is always non-degenerate. Moreover, because cr^r*^ 
commutes with ct'^t^, we have masses 777 — A and 777 -I- A 
for the Dirac fermions. Explicitly, the spectrum is 
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Note that one of the masses changes its sign at m = A, 
thus there is a QPT. Except at the QPT point, we have 
a fully gapped spectrum. Upon changing the sign of the 
mass, the system could change its chirality. To see that it 
does, we do a unitary transform to get a block-diagonal 
form of Q: 
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Erom this representation, it is clear that we have a QPT 
at 777 = A as there is an abrupt change in the winding 
numbers at 777 = A. Each block has the winding number 




FIG. 2: Illustration of the dependence of band structures on 
the parameters m > and A > 0. Here, we assumed A 7^ 0. 
There are always 4 bands, and the bands (blue) with the 
gap m -|- A never touches each other. The bands (red) with 
the gap A — m touch each other at m = A, signaling the 
quantum phase transition. A. The band structure for the 
case A < m. The system is gapped. B. When A = m, 
the spectrum becomes gapless and there is a change in the 
winding number. C. For A > m, the system is gapped again. 



sgn(m + A)i and sgn(777 — A)^. For |A| > \m\, we have 
two cones of different chirality and thus we are in QSH. 
However, if we pass through | A| < 1 777], we have two cones 
of the same chirality and thus we are in QAH (See Fig[2|. 



Note that our QSH phase does not support the helical 
edge state at the boundary of the sample (however, see 
the appendix) though the quasiparticles in the bulk low- 
energy theory can be considered in the standard QSH or 
two-dimensional topological band insulator phase. When 
the higher order kinetic energy correction to the Dirac 
Hamiltonian is included (e.g the hexagonal warping term 
^ ^{k\ -\- k^)a^), the winding numbers for QSH phase 
are changed from ±1/2 to zero for the both Dirac sectors. 
However, the low energy theory of our QSH phase is the 
same as that of the standard QSH physics. 

In this section, we have seen from the band structure 
of the Hamiltonian ([T]) that there is a QPT between QSH 
and QAH phases. The QPT is driven by the competition 
of magnetic mass 777 (which prefers the QAH phase) and 
exciton mass A (which prefers the QSH phase). At the 
transition, the mass gap vanishes and there is a change 
of the band structure's chirality. Except at the transi- 
tion point, we always have a gap for the single-particle 
spectrum. The goal of the next section is to obtain the 
key universal properties of the two phases in an abstract 
approach that is less dependent on microscopic details. 
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III. EFFECTIVE THEORY FOR VORTICES 
AND FLUCTUATIONS 

Now that we have identified the chirafity and chern 
numbers, we integrate out the fermions to obtain efltective 
descriptions for the fluctuations and vortices of the order 
parameters and gauge fields. First, we identify which 
fields and vortices might appear in the effective theory. 
To do this, we include all possible interactions and write 
them as gauge fields acting on the Dirac fermions. We 
also include Lfci„(A), is the kinetic term for A. 

a 

+ E(/o^Aa"/i + M+ife.n(A) (10) 

Here, A" is the electromagnetic gauge field of the layer 
a — 1,2. From the gauge fields and , we can de- 
fine A = A^ + A"^ which couples to the total electromag- 
netic charge and /3 = {A^ — A'^)/2 which couples to axial 
charge. Before getting into the details of the effective field 
theory, we study the nature of the vortex in the conden- 
sate. Especially, whether the vortex is accompanied with 
the gauge flux or not is crucial for the properties of the 
vortex. So, we begin with the kinetic term Lfci„(A) for 
the exciton condensate. Due to the gauge symmetry of 
the exciton order parameter A =< f^^f^ >, the axial 
gauge couples to A minimally 
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which is XY - model (of the phase stiffness Ps — |Ap/C/) 
coupled to the gauge field /3^. Now, we imagine the vor- 
tex configuration of which can be generated by the 
solenoids placed near to the exciton condensate (Figjsj). 
The simil ar exp erimental scheme has been considered in 
the paper^i^mi. Then, the standard dual transformation 
allows us to write — 2/3^; — for the vortex config- 
uration of the exciton condensate. Hence, § P = ttZ for 
the vortex in the condensate. This duality can be readily 
seen if we do the particle-hole transformation of the exci- 
ton order parameter A =< /^Vct > fo^' the layer 2 only. 
Upon the particle-hole transformation of the layer 2, we 
see the exciton order A becomes the superconducting 
order parameter (A — )•< fl^ fl^ >). Thus, the exciton 
condensate is a 'superconductor' in the axial channel and 
an insulator in the total charge channel as noted befor^^^. 
Here, if we have an axial vortex configuration for (3^, we 
have the counterflow (or an axial current) Jcf ^ Psdfj.4' 
encircling the vortex to screen the phase gradient due to 
(or vice versa). We estimate the screening length by 



||a^0-2/3^p-f— (9/3)2 + —(9^)2 (12) 



where the last two terms are the kinetic energy for the 
gauge /3 and A which is obtained from the maxwell term 
for Aj^ (= + A^/2) and Al (= - A^/2). Hence, we 



identify an axial screening length As ~ l/^/p^. Note 
that there's no screening length for total electromagnetic 
gauge field which reveals that we have an insulator 
in total charge channel. In principle, there should be an- 
other term ~ (/J^/ed)^ in Eq.(12) when the film height 
d gets thin enough. However, we ignore this term upon 
assuming that the film is not too thin d > 0{1). Though 
the energy of the vortex in the exciton superfluid could 
be reduced from infinity (~ log L of the system size L) to 
finite value by screening, it might be costly for the gauge 
field Pfj, for the vortex configuration. We will compare the 
energy cost for the gauge configuration of the screened 
vortex with the energy of the unscreened vortex and see 
when the vortex will tend to be screened by /3^. To do 
so, we roughly estimate the energy cost for the magnetic 
field configuration spontaneously driven (without the ex- 
ternally imposed magnetic fields in the condensate) by 
the screened vortex ^ /3 = tt. We imagine that the mag- 
netic field is dragged from infinity to the center of the 
circular exciton condensate of radius L (with height d) 
and pulled out through the vortex positioned at the cen- 
ter of the condensate (See figures B. and C. of Fig|3|. 
Magnetic field configuration B{r) = B{r)er is 



B{r) 







27rr • d 



(13) 



with the unit vorticity $o- The electromagnetic en- 
ergy for this configuration / ^XeB^ grows as Xe^Q^ogL. 
Hence, if we compare the energy ^ ps^plogi of the 
unscreened vortex and the energy ~ Xe^o^ogL of the 
screened vortex, we conclude to have screened vortices if 
the magnetic permiability of the exciton condensate Xe 
is much smaller than the phase stiffness ps ■ When this is 
not the case, there could be more interesting possibilities 
such as the irrational charge and statistics for vortices. 
We refer the reader for this discussion to the papei!^ and 
the references there-in. On the other hand, if the mag- 
netic field is supplied by the external mean such as the 
experimental set-up in our scheme (FigjsJ) or Fig.l of the 
paper Jl^Ell^ it's much easier for the condensate to have 
the screened vortices. The vortex in the condensate can 
save logarithmically divergent energy cost (^ log L) by 
slightly tilting and trapping the magnetic fluxes which 
should take finite energy, and this concludes that the 
vortex in our experimental scheme will be always accom- 
panied with the gauge flux §13 — ttZ (FigjsJ). 

Another comment is in order. Even though the exci- 
ton order parameter A in the mean field theory is of the 
same form to the tunneling term between layers, tun- 
neling term introduces a new term to the effective field 
theory Li„„ ~ — tcos(0) as it tends to pin down the U(l) 
phase of the exciton order parameter [t is proportional to 
the tunneling strength between layers) . This term breaks 
U(l) symmetry and induces 'sine-Gordon' theory coupled 
to the vortex matter which can be ignored as i — )■ in 
our perspective^. 

Now that the allowed gauge couplings are identified, 
the fermions are integrated out to obtain the effective 
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FIG. 3: Illustration of the generating axial gauge fluxes. Gen- 
eral idea is to make magnetic field configuration which has 
J (BI — B1) nonzero. A. Experimental scheme similar to that 
of the papei^. We place magnets or solenoids at the bound- 
ary of the sample. B. Experimental scheme similar to that of 
the papei'^. We place solenoids on the top and below of the 
condensate. Without the solenoids, the vortex should create 
the magnetic flux itself to be screened. To estimate the en- 
ergy of the vortex roughly (see the text), we consider that the 
flux enters at the vortex and leaves at the boundary of the 
sample. C. Top view of the scheme B. The exciton condensate 
has the radius L, and the solenoids are placed at the center 
of the condensate. 



theory which is the main result of this section: 



L = sgn{m + A) — {A + l3)^d,{A + I3)x 

+ sgn{m - A)^(A - P)^dM " P)x + 0{{dAf) 
on 



(14) 



(with e — 1). Here, we ignore the Maxwell term as 
the Chern-Simons (CS) term dominates the low energy 
physics. The appearance of the doubled CS theory is due 
to the two massive Dirac cones in Eq. ([7| where we gauge 
one Dirac cone of mass {m + A) with A^ - ^ 
other Dirac cone of mass (m — A) with Af^ 
tegrate out those gapped fermions to obtain the gauge 
theory of and A^. This phenomenological approach 
reproduces the essential feature of Eq.(14) which is con- 
sistent with the previous studji^l. We will also justify this 
effective theory Eq. ( 14 ) in the later part of this section. 



and the 
/3fj. and in- 



However, readers interested in a more detailed derivation 
of Eq.p^ and Eq. 



!|18| from the Dirac hamiltonian are 
referred to the papei^ where this field theory is rigorously 
considered. In this field theory, we observe the coefficient 
of the CS term for A changes aX m ~ ±A. This signals 
the QPT at m = ±A. This simple theory captures nearly 
every essential property of the vortices and fluctuations 
of the system; the rich behavior of the vortices includes 
localized charges. 



For the QSH phase(|A| > |m|)), we have (up to the 
sign factor infront of the mutual CS theory for and 
A,) 



Zn 



(15) 



From this effective mutual CS (or equivalently BF) the- 
ory for and A^, we can directly read off charge and 
statistics information for the exciton vortices. Here, it's 
interesting to notice that BF theory emerges for QSH 
phase which reflects the time reversal symmetric nature 
of the phase. We now study the charge and statistics 
of the vortices. First, we see that the absence of the 
topological term for /3 requires that the vortices have 
no statistics (statistical angle 9 — 0). However, the 
EM charge localized at the vortices is non-trivial in that 

9 = / M = / 2¥^/5 ^ h§l^ ^ ±1/2, that is, we have 
a half of electron charge localized at the vortex, which 
is consistent with the previous consideration on the vor- 
tices of TEC ( A = 0) . In comparison with the papeil^ and 



our result Eq.(15), we note that we miss the coefficient 
sgn[^s) for BF term which matters for the sign of the 
fractional charge at the vortex but this sign factor is not 
important for our thin film topological insulator prob- 
lem. The electric charge of the fermionic zero mode at 
TEC is known to be defined modular integei^ in connec- 
tion with the axionic 6*- vacuum i.e., the localized charge 
at the vortex is — e(l/2 -|- n) for n € Z, and thus the 
sign of the fractional charge (or the sign of BF term in 
Eq.(15)) shouldn't be taken seriously. As was noticed in 
the earlier study of the fractional charge, the ±e/2 charge 
is associated with the zero mode solution for the vortex. 
Therefore, we look for the solution of H^* = such that 

H = 
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ujpE^. Then, the single phase method provides 
Uo ~ e~*^ and vi,V2 independent of the angular 



(16) 

which is the Hamiltonian for the exciton vortex with vor- 
ticity —1. We solve this problem by the following ansatz 
^ ^ {ui , Vi, V2, U2)'^ with the constraint u\ — —U2 and 

Ui ^ u>2 _ ^ _ 

variable With this in mind, we use the following 

ansat^ 

ui(r,6l) = /(r)e"/4e-■/"■'^('■')'''■'e-'^ 

„i(r,6l) = (7(r)e-*^/'*e--'"^^('^')'''-' (17) 

Upon substituting this ansatz, it is straightforward to 
solve the resulting differential equation and the so- 
lution is u\ ^ exp{— A[r')dr')Io{mr) and vi ^ 
exp{— A{r')dr')Ii{mr) where /„ is the nth modified 
Bessel function. This wavefunction is convergent if and 
only if we are in QSH phase |m| < |A|. 
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The statistical phase and locahzed charge of the vortex 
can be understood by the following adiabatic argument. 
In the QSH phase, we have two gapped Dirac fermions 
and they carry half-filled quantum Hall effect with the op- 
posite chirality. When we thread the axial flux <f /3 ^ tt 
adiabatically to the system which is identical to an ex- 
citon vortex of the circulation 27r, each quantum Hall 
state collects charge e/4. Thus, the localized charge of 
the vortex is e/2 (= e/4 + e/4). When we do the pair- 
wise exchange of the vortices, the statistical phase accu- 
mulated is (= 7r/8 — 7r/8) because the statistical phase 
7r/8 from one quantum hall state is precisely cancelled by 
the phase — tt/S from the other. Hence, the vortex in the 
QSH phase has no statistical angle: 6 = 0. 

In terms of the external EM gauge A, we have the 
response ^ 0{dA)'^ which is the Maxwell kinetic term 
and non-topological. However, the topological property 
shows up in principle if we consider a particular type of 
the domain wall for the exciton mass A. (Appendix) 

For the QAH phase(|A| < |m|), we have 

L = ^e'-^^A^d^Ax + ^e^'-'/S^^.^A + OidAf (18) 
in An 

From this effective theory, we now notice that we are 
in the quantized QAH from the leading CS term for the 
response of A. This tells us that we will have a circulating 
chiral edge state supporting the Hall conductance e'^/h 
at the boundary of the sample. 

On the other hand, we have no crossing term between 
A and /3 implying that the vortices of the exciton order 
will not carry any EM charge. On the other hand, the 
exciton vortices is ^anyonic\ It looks like the coefficient 
of CS term implies that the vortices are fermionic, but 
this is not the case as we have §13 = ±7rZ (In the usual 
CS theory L = -^ada, § a = 2tt defines the unit vortex). 
With this consideration, we now have fractional statis- 
tics between vortices with the statistical angle 9 = tt/A 
(half of semionic statistics) . Note that the vortex in QAH 
phase does not carry a zero mode solution for the equa- 
tion (16 1. However, we can still talk about the exciton 



vortex as it is quantized to have circulation 27rZ and 
well-defined excitation. This anyonic vortex emergent in 
QAH phase can be considered as an example of anyons 
from 'weakly' interacting systems and this is similar to 
the previous study on integer quantum hall effect with 
the fully filled lowest landau level adjacent to the type- 
// superconducting filnP. 

The behavior of the vortex can be understood in terms 
of the adiabatic argument as before. In the QAH phase, 
the two sectors of gapped Dirac fermions have the same 
chirality. Thus, the axial gauge flux <f /3 = tt associated 
with an exciton vortex will collect e/4 in one sector and 
—e/4 in the other. This gives the total charge (= e/4 — 
e/4). If we exchange a pair of vortices, we now have the 
statistical phase 7r/4 (= tt/S + tt/8) by adding up the 
phases accumulated from the two sectors. 

We now take a different point of view to present the 
justification of our theory Eq.(14|. We base on the adi- 



abatic argument and physical understanding of the frac- 
tional charge and statistics of vortices. First, we expect 
that there could be (at most) two CS terms for the effec- 
tive field theory of /3p (axial gauge) and (electromag- 
netic gauge) due to the two massive Dirac cones 



^eff - 



2n 



e^^'^^A^d,l3x+^e^''^A^d,Ax+^e''''^f3^d,/3^ 



(19) 

and the effective field theory should behave well because 
of the gap except the critical point, i.e., the coefficients 
(Ci,C2,C3) in Eq.(19l change only at the critical point 
TO = A. So our basic strategy is to extract the coefR- 
cients for QSH and QAH phases. As T-symmetry has the 
important role, we study the transformation of and 
Afj_ under T-symmetry operation: (/3o,/?i) ^ {—(^OiPi) 
and {Ao,Ai) — >■ {Ao,—Ai). Thus, the first BF term in 
Eq.(|l9|) between and A^, is T- symmetric, and the 
other two CS terms are T-breaking. 

For QSH phase, the phase should respect T-symmetry 
in the effective field theory Eq.(19) though we break T- 



symmetry microscopically by inclusion of magnetization 
mass ^ TOcr^. We can extract the coefficient (Ci, C2, C3) 
for Eq.(19) at to = with the fixed A because the QSH 
phase with finite \m\ < |A| is adiabatically connected to 
the system of to = without closing the gap. This con- 
sideration gives C2 = C3 = for the entire range of QSH 
phase. And we supplement the field theory with the in- 
formation that the vortex § 13^ = it carries the fractional 
charge ±e/2 with finite A at m — 0^. This translates to 
the field theory as Ci = 1. Hence, we have the effective 
field theory for QSH phase L^f f = ^ f3dA (with the im- 
plicit antisymmetrization of indices for fields) and this is 
precisely the same as the field theory we obtained before 
Eq.pl). 



On the other hand, QAH phase has the total Chern 
number which fixes C2 = 1. For Ci and C3, we utilize 
the adiabatic argument with A — while keeping to 
finite. At this limit, the axial gauge flux § (5 = n as- 
sociated with an exciton vortex will collect e/4 in one 
layer and —e/4 in the other. This gives the total charge 
for the vortex which translates Ci = for the effec- 
tive field theory Eq.(19). When we do a pair- wise ex- 
change of vortices § /3 = tt, we now have the statistical 
phase 7r/4 (= tt/8 + tt/8) by adding up the phases accu- 
mulated from the two layers. This fixes C3 = 1. As a 
whole, we obtain the effective field theory for QAH phase 
L = j:^AdA+ -^(38 f3 (with the implicit antisymmetriza- 
tion of indices for fields) which is the same as in Eq.( 18 1. 
Finally, it's straightforward to rewrite these field theories 
in both phase into a sinlge theory of the form Eq.( 14). 

In the papei^i^, it's shown that the excitonic vortex 
may carry the irrational statistics due to the irrational 
axial charge SQy = — bound to the vortex. (Note 
that in the bilayer system, the fluctuation in the axial 
charge could be finite unlike the total charge which is 
suppressed by the Coulomb interaction between layers) 
With this additional axial charge and the axial tt flux of 
the vortex, the statistical angle is given as ttSQ^ which 
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could be continuously tunecP^. In our case, this might 
not be the case. First of all, there's no term breaking 
the symmetry between layers (such as the voltage drops 
/is), we don't have any irrational statistics as the axial 
charge is SQv = due to the symmetry between up and 
down layers. We also can deduce the form of the effective 
field theory from the adiabatic argument when the volt- 
age drop fis between layers is included in consideration. 
When the voltage drop fis is slowly turned on from in 
the system, we evolve the mass A in Eq.( 15 1 and Eq.( 18 ) 



into + without closing gaps (the total masses for 
the two Dirac cones are m ± -y/ A^ + /i^). So we predict 
at least in our effective field theory and the adiabatic 
argument that the statistics and charge of vortices are 
fractionalized and quantized as in the case without fig . 

In the QAH phase, the system acquires a topological 
QAH response to the external EM gauge. This has a 
direct implication on the Faraday angle Op which is re- 
lated to Hall conductance a^y of the systenPlMini, jn 
QAH phase, we expect to have a strong Faraday effect as 
T-symmetry is broken. On the other hand, we wouldn't 
have a strong Faraday effect in QSH phase as we do 
not break T-symmetry effectively though we break T- 
symmetry microscopically. This optical response could 
be used in experiment to distinguish the two phases. In 
a thin film with interlayer coherence (exciton order) , the 
electron tends to be ambiguous on its layer index. Thus, 
we can think of the thin film as like a single layer with 
two Dirac cones. Then, as pointed out from the previous 
section, we have two Dirac cones with masses m-\- A and 
TO — A, and the two Dirac cones contribute to total a^y 
equally. We specialize for low frequency uj << Ec {E^ is 
the cut-off for the energy, and typically we can take Ec 
as the bulk gapSl) and the chemical potential lying in the 
gap. Then, at the leading order 0{uj / Ec,m/ E^ A./ Ec) 



0, ax 



(20) 



where (Jxy,± = sgn(TO ± A) x ^(1 - ) (Here, a 

in the coefficient is the fine constant, i.e a — 1/137). 
For convenience, we take the Zeeman mass to and exci- 
ton order parameter A positive. In the QAH phase, we 



have a 



xy - 

tan 



' aj^Tx with the limit m/Ec — 
/ ^" / ) 10^3 rad. 



0. This gives 
which is sim- 



" V<^7a''-i-\/'^/'^ 

ply the double of the previously studied on the single- 
layer Dirac cone with T- breaking. In the QSH phase, 
where we preserve T-symmetry effectively, we have to- 
tally different behavior in that axy = x ^ — in the 
limit ra/Ec -> 0. Thus, there is no significant Op even 
though T-symmetry is broken at a microscopic level. The 
small T-breaking shows up only in the order of 0{m/Ec)', 
9p ^ axO{m/Ec). If we plug to ^ 10 meV and Ec ^ 0.3 
eV for Bi2Se3, then we have Op ^ a x 0{m/Ec) < lO""* 
rad, much smaller than Op for QAH phase. 



IV. UNIFORM ORBITAL MAGNETIC FIELDS 
IN A THIN FILM 

We have seen that the interplay between Zeeman and 
exciton masses induces interesting physics on the vortices 
and in the electromagnetic response. However, there is 
another natural way to break the T- symmetry: uniform 
magnetic fields along z. We have Zeeman interaction due 
to magnetic field B via ^ gcr^B but we also have Landau 
Levels (LL). We will see that the zeroth LL will deter- 
mine the physics and we can obtain the same effective 
theory as before. The appearance of the same effective 
theory can be traced back to the CS effective theory for 
QHE. Now, the coefficient of the effective theory is de- 
cided by the filling of the LLs, rather than the winding 
numbers of the band structures, and the phase transi- 
tion at TO = A is replaced by the quantum hall phase 
transition where the filling of LLs are suddenly changed. 
We begin with LLs of Dirac fermions which is similar to 
graphene. But the crucial difference here is the Zeeman 
coupling and degeneracy. Under the uniform magnetic 
field B, electrons form the LLs with index N G Z 



E{N,r) = sgn{N)\vp\y/C\N\ + (to + sgn{r)A)^ 

r^±l,N^O (21) 



and C 



ehB^ 



The zeroth LL is sensitive to the compe- 



tition between to and A in that (See also Fig. [4| 

E{N = 0,r) = m + sgn{r)A, r = ±1 (22) 

We begin with the case where the chemical potential lies 
at Ep = as before. Then, we see that every LL with 
negative (positive) index A'^ is always filled (empty) in- 
dependent of the parameters B, to and A. However, the 
filling of the zeroth LL is dependent of the parameters 
(See FigQ. For to > A, we have that both of zeroth 
LLs are empty (as they are above the chemical potential 
fi = 0). Thus, we have QHE with the quantized hall con- 
ductance ^ (Note that there is an offset by 1/2 of the 
QHE for Dirac fermions). On the other hand, if we have 
TO < A, one LL is filled and the other is empty so we 
have a spin-Hall-like response. For the QHE, CS theory 
is the effective field theory 



L, 



QHE 



AttK 



(23) 



where K is the inverse of the filling (i.e K — 1/v). Hence, 
we conclude the effective theory for the system as K = 
±2 (The sign of K is determined by to and A) 

Lqhe = sgn{m - A)-^Aif_,d^Aix 



Sit 



+ sgn{m + A)—A2^d,A2x + 0{{dAY) 



(24) 



Stt 



This is exactly the same effective field theory before ( |14[ ) , 
and what all we need to do is to assign the correct cou- 
pling of axial gauge /3 for each of the zeroth LLs. The two 
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A. 



B. 



m 



With this ansatz, we obtain four independent hnear equa- 
tions of f±{r) and g±{r) 



d_ 
dr 
d_ 
dr 
d_ 
dr 

or 



mf+ + A.g_ + ( 
mg+ + A/_ + ( 
m.9_ + A/+ + ( 



r 2 

r 2 
Br, 

-^)9- = 



(26) 



We can try to solve this system of differential equations 
by reducing into the following two equations via canceling 
/_ and in the above equations. 



FIG. 4: Illustration of the Landau Levels depending on the 
parameter m > and A > 0. Here, we assumed A 7^ 0. 
We always keep the chemical potential ^ at iS = 0, and the 
0th LL passes through the fermi energy i? = at m = A. 
The levels of LL index n with n 7^ never change their filling 
as their energies are always positive (if n > 0) or negative 
(if n < 0). These levels (blue) do not play an important 
role. On the other hand, the two 0th LLs drive the quantum 
phase transition at m = A because the sign of their energies 
are dependent on m and A. The filling of the lower 0th LL 
(red) in the figure is abruptly changed at m = A as it passes 
through /i. A. the LLs of the case A < m. The lower 0th 
LL is fully- filled B. when A = m, the 0th LL experiences 
change in its LL filling. C. Passing A > m, the 0th LL is 
fully emptied 



LLs are from the two Dirac cones of the opposite axial 
charge ±1 and the same EM charge -1-1, i.e., Ai = A + /3 
and A2 = A — /3. Upon plugging this in, we restore the 
same effective field theory Eq( 14 1 as in the previous sec- 
tion. Thus, we have the same charges and statistics if 



we include the vortex fields in the effective theory. We 
see immediately that much of the physics studied in the 
previous section applies in this case. 

However, the zero energy solution for the magnetic 
field case is not trivial to generalize as the uniform mag- 
netic field cannot be gauged away. Furthermore, this 
unifrom orbital magnetic field destroys the particle-hole 
symmetry of the Hamiltonian, and thus we cannot use 
the same ansatz Eq.(17) with ~ (ui, ui, ti2, 1*2)^ with 
the constr aint u* = —U2 and = V2- This form of 
the ansatJi^E^ relics on the anticommutation relation 
{C,H} = where C = iTVayK and K is the com- 
plex conjugation. With the orbital magnetic field, we no 
longer have the anticommutation relation {C, H} = 0. 
This {C, i/} 7^ manifests in the LL spectrum (See 
Figj4|. So we instead make an 'educated guess' for the 
ansatz (m+, — w_)^ and 



(A/2 + 92 _ 
(M^ 



9 D d2^2 



B B^r 



2^2 



m{2d 



= d and AP = 



where we abbreviate 
above differential equation Eq. ( 27 



- + Br)g+ = 

- - Br),U = 

(27) 

-A2. The 



is difficult to solve and 
not clear if they admit analytic so utions without approx- 
imations. Rather than directly attempting to solve the 
system of differential equations, we look for a few solvable 
limits with the analytic solutions in terms of the conflu- 
ent hypergeometric functions. The most convenient and 
important limit is when the zeeman coupling vanishes 
(to —7' 0) which effectively decouples f± and g±. In this 
limit, the solution for Eq.(27) is /± = and 



9+ir) 



M 2B 



1 Br^, 
4' 



(28) 



where VF( ^ - ^ ; - 1 ; ^ ) is the Whittaker functioiP, or 
equivalently we can represent it as the parabolic cylinder 
function g+{r) ^ Dp{^^-) with P — —jg which appears 
in the Kekule vortex solution in the graphen^^fll under the 
uniform orbital magnetic field. This limit of the solution 
corresponds to the case |A| >> |to|, i.e., the negligible 
zeeman coupling. 

The other important limit where we can obtain an an- 
alytic solution for Eq. (26) is when r ^ 00. In this limit, 
we ignore the potential terms of the order 0(l/r) for 
f± and g± in Eq. ( 26 1 . From this approximation, we try 



to see if there's a convergent solution for the system of 
differential equations. In the limit r — > 00, we obtain 



9+{r) 



(A - m'f 



1 Br^ 



2B 

(AH 



) 



2B 



2 

1 Br^ 
4' ~Y~ 



)} (29) 



u±{r,9) = /±(r)e±*'^/''e--'"^'^(''')'^'^'e~'^ 



v±ir,0) = g±{r)e 



=Fi7r/4g- A{r')dr' 



(25) 



with the initial condition dependent coefficients Ci and 
C2. It's now obvious that the solution g+{r) in Eq.(29) 
reduces to g+{r) from Eq.(28) when to/|A| — 0. As 
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9+K 



„2M 



exp(— i3r^/4)(l + 0(l/r^)) as r — ^ cxo where 
M = (to ± A)^, the ignored potential term ~ 0(l/r) 
in Eq. ( 26 1 wiU not generate more singular terms than 

^2M 



exp(— _Br^/4) in power series expansion near ; 
Note that the exciton vortex in the uniform orbital mag- 
netic case carries the same charge and statistics as the 
case without the orbital field due to the same form of the 
effective topological field theory despite of the different 
forms for the zero mode solutions. 



e* between two topological insulators (See FigjljB), and 
which will change — ^ e^/e*f/ which could be as large 
as 3 if e* ^ 10. This might open a way to achieve a 
dramatic increase in the superfluid transition tempera- 
ture and realize the surprising physics of the topological 
exciton condensate. 
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V. CONCLUSION 

In summary, we have studied a thin film of topologi- 
cal insulator with both T-breaking Zeeman mass and T- 
symmetric excitonic mass. The two masses compete with 
each other and result in two topologically distinct phases 
for the elementary excitations: quantum anomalous Hall 
(QAH) and quantum spin Hall (QSH). We studied the ef- 
fective theory for the electromagnctism and exciton order 
parameter vortices by integrating out the fermions, and 
there can be other topological properties such as a helical 
metal at a particular kind of domain wall of the exciton 
order parameter. We also obtained explicit wave func- 
tions for the fermion zero mode at the vortices. There is 
one zero mode for QSH and no mode in the QAH regime, 
and found the zero mode solution under a uniform mag- 
netic field in some limits. 

Before finishing this paper, we would like to empha- 
size some relevant facts for practical observation of ex- 
citon condensation in a topological insulator (TI) thin 
film. In principle, TI might host a higher transition tem- 
perature T* for the exciton superfiuid-insulator transi- 
tion than graphene, resulting from the decreased num- 
ber of fermion flavors. In graphene, there has been in- 
tense theoretical study of the exciton condensates, and 
the estimated transition temperature for graphene ranges 
from milliKelvirP2! to room temperature?^]. An impor- 
tant reduction comes from the large number of fermion 
species in the screening^^ [N = 8 for bilayer graphene) 
which induces the factor T* ~ e~^^Ef where Ef is the 
fermi energy of the system even though the bare in- 
teraction energy scale (without screening) is not small 
Ts = e^/evf ~ 1 in the graphene. However, the thin-film 
topological insulator contains only two species N = 2 
because we have only one Dirac fermion per layer. But 
we have the disadvantageous situation for as in many 
current materials the dielectric e for TIs are quite large 
(ranging from 30 to 80) though f/ ~ 10^to/s which is 
nearly half of graphene case. (An additional problem of 
early TI materials, that they were not in practice very 
insulating, seems to have been overcome.) These two 
factors contribute to reduced values ^ 1/3—1/8. How- 
ever, we can look for higher and thus higher transition 
temperature by noting that the main drawback comes 
from the large e which is material-dependent; a possible 
solution is to sandwich a thin normal insulator of small 
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APPENDIX 

We show that if the exciton mass has a "real" domain 
wall (i.e., passes through zero with some fixed phase), 
there is a topological helical state. Without loss of gen- 
erality, let the exciton mass A(y) Aq > for y — > — oo, 
and A — — Aq < for ?/ — )• oo with A(?;) — >• as y 0. 
Then we have a helical edge state localized aX y = 0. In 
terms of the four-spinor, the helical states are 5'+(fc) ^ 
exp{ikx) exp(-| A{y)dy\){l, 1, 1,-1)^ with E{k) = k 
and *_(fc) ^ exp(zfca;)exp(-|jyA(y)dy|)(l, -1,1,1)^ 
with E{k) = — fc and two fermionic states 4'+(fc) and 
^_(— fc) are the Kramer pair. These ID state is protected 
when T- symmetry at the domain wall is conserved as for 
spin Hall edge states. Thus, this helical state localized at 
the domain wall reflects the underlying topological states. 
This spin Hall type physics can be best understood if we 
look at the matrix form of the Hamiltonian H = ^f^H^ 
with the proper unitary transformation from the original 
bases. 



H 



- A 


P+ 





- 


P- 


-A 














-A 


P+ 








P- 


A _ 



(30) 



As the helical ID metal is protected only if T-symmetry 
is conserved at the domain wall, we set to = near the 
domain wall. Now, it's clear that when A has the domain 
wall at J/ = 0, the upper two spinor has the chiral mode 
propagating in the positive x, and the lower two spinor 
has the chiral mode propagating in the negative x. Fur- 
ther consideration shows that two modes are a Kramer 
pair and are protected from generating a gap only when 
T-symmetry is conserved. Note that a general domain 
wall where the phase difference in the (generally com- 
plex) exciton mass is not tt need not have this bound 
state. 
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